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The  U.S.  Army  Corps  of  Engineers,  Seattle  District  has  completed  a  feasibility  study  and  determined  barrier  island 
restoration  to  be  the  most  appropriate  long-term  coastal  flood  and  storm  damage  reduction  measure  for  the  Shoalwater 
Indian  Reservation.  Over  the  last  century,  Cape  Shoalwater  has  receded  more  than  2.8  miles.  By  1990,  the 
Shoalwater  Reservation’s  only  remaining  protection  from  storm  wave  attack  was  a  series  of  barrier  islands  fronting 
Tokeland  Peninsula.  Extreme  water  levels  coincident  with  strong  winter  storms  have  historically  inundated  this  low 
lying  topography  and  are  responsible  for  the  erosion  and  overwash  of  the  protective  barrier  island  known  as 
Graveyard  Spit.  Here  a  simple  risk  assessment  tool  is  presented  for  identifying  flood  risk  to  the  Shoalwater 
Reservation  infrastructure.  Statistical  analysis  of  extreme  water  levels  and  numerical  modeling  is  utilized  to 
determine  the  extent  of  inundation.  From  the  analysis  it  was  determined  54%  of  the  inventoried  infrastructure  is  at 
risk  during  a  storm  event  equivalent  to  the  observed  event  on  March  3,  1999.  With  the  barrier  island  restoration  it  was 
found  that  this  risk  is  reduced  to  7%. 
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INTRODUCTION 

The  Shoalwater  Bay  Indian  Reservation  is  located  on  the  northern  shore  of  Willapa  Bay,  a  large 
natural  estuary  located  on  the  southwest  coast  of  Washington  State,  USA.  With  a  diurnal  tidal  prism 
of  more  than  2.5  x  1010  ft3  (Jarrett  1976),  Willapa  Bay  is  one  of  the  largest  inlets  on  the  coast  of  the 
continental  United  States.  The  northern  shoreline  of  Willapa  Bay  has  changed  substantially  since  the 
Shoalwater  Reservation  was  established  in  1866  (Terich  and  Levensellar  1986).  Over  the  last  century, 
portions  of  the  Cape  Shoalwater  shoreline  have  receded  more  than  2.8  miles.  By  1990  the  Shoalwater 
Reservation  only  remaining  protection  from  storm  wave  attack  was  a  series  of  barrier  islands  fronting 
Tokeland  Peninsula  (see  Figure  1).  Tidal  currents  and  storm  waves  have  severely  eroded  the 
protective  barrier  islands  known  as  Graveyard  Spit.  A  series  of  breaches  beginning  in  1 994  exposed 
the  Reservation  to  increasing  levels  of  flooding  during  extreme  water  levels.  Erosion  from  storm 
overwash  has  filled  North  Cove  with  sand  and  significantly  degraded  and  altered  the  habitat  of  the 
once  protected  embayment.  The  increased  wave  energy  through  North  Cove  has  led  to  an  increase  in 
the  severity  and  frequency  of  flooding  and  erosion  of  Reservation  uplands  during  storm  events.  In 
2000,  the  Assistant  Secretary  of  the  Army  for  Civil  Works  authorized  the  U.S.  Army  Corps  of 
Engineers  to  construct  and  maintain  a  project  at  Federal  expense.  The  feasibility  of  hard  and  soft 
engineering  alternatives  were  evaluated  by  a  multi  agency  technical  team  and  the  analysis  determined 
barrier  island  dune  restoration  with  periodic  nourishment  was  the  most  appropriate  long  term  solution 
to  mitigate  future  erosion  and  flood  damage  risk  (US ACE  2009a). 

SITE  CONDITIONS 

Geology 

Graveyard  Spit  exists  as  a  thin  and  fragmented  landform  that  is  anchored  and  aligned  by  consolidated 
and  erosion-resistant  Pleistocene  substrate  in  the  vicinity  of  the  SR  105  emergency  stabilization  groin 
(Figure  1).  Extending  to  the  east  of  Graveyard  Spit  is  a  series  of  segmented  sand  islands.  In  contrast 
to  historical  conditions,  this  fragile  line  of  barrier  beaches  no  longer  appears  to  receive  sand  supply 
from  the  eroding  beach  plain  to  the  northwest.  The  lack  of  sand  supply  indicates  that  this  landform 
will  remain  of  low  relief,  compromising  its  historical  function  as  a  flood  barrier  for  the  Tokeland 
Peninsula. 

Shoreline  retreat  along  this  northwest  comer  of  North  Cove  has  slowed  substantially  relative  to 
historical  rates  of  change,  but  the  present  condition  and  orientation  of  the  spit  suggest  that  it  will 
continue  to  pivot  towards  the  north-northeast  from  its  hinge  point  at  the  base  of  the  Pleistocene  terrace. 
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Thus,  the  present  condition  of  the  spit  is  locally  controlled  by  the  geological  framework  of  the  region 
(USGS  2004) 

Tidal  Currents,  Water  Levels,  and  Wave  Climate 

Tidal  currents  are  strong  in  the  north  entrance  channel  which  is  presently  the  main  conveyance  at 
the  mouth.  The  primary  entrance  channel  at  the  mouth  has  migrated  from  east-west  to  a  northwest 
orientation  since  1941  (Kraus  2000).  Peak  spring  ebb  currents  through  the  north  entrance  channel  are 
approximately  3.3  knots  while  peak  flood  currents  are  approximately  1.6  knots. 

The  mean  tidal  range  at  the  entrance  to  Willapa  Bay  is  6.8  ft,  and  the  peak  spring  tide  elevation  is 
approximately  1 1  ft  above  mean  lower  low  water  (MLLW)  as  published  by  the  National  Ocean  Service 
(NOS)  of  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  at  the  Toke  Point, 
Washington  Station  9440910.  In  winter  months  spring  tides  are  coincident  with  large  storm  events  and 
are  capable  of  producing  water  levels  exceeding  5  ft  the  predicted  astronomical  tide.  Table  1  lists  the 
maximum  observed  water  levels  at  the  Toke  Point  Station  since  1973.  Storm  durations  during  such 
events  range  from  25  to  124  hours  as  shown  in  Figure  2. 

Field  measurements  from  a  nearshore  instrument  deployed  offshore  of  Graveyard  Spit  in 
December  2002  measured  peak  significant  wave  heights  of  approximately  3-4  feet  during  high  tides 
(US ACE  2009b).  These  field  data  confirmed  wave  energy  is  significantly  dissipated  by  the  offshore 
shoals  at  the  entrance  and  the  wave  energy  propagating  over  Graveyard  Spit  is  tidally  modulated. 


Figure  1.  Shoalwater  Bay  -  Barrier  Island  Restoration  Project,  Willapa  Bay,  Washington  USA 


Table  1.  Toke  Point  Highest  Tides,  1973  through  2010 

TOTAL  WATER 
ELEVATION 
(FEET,  MLLW) 

MAX  SURGE 
ELEVATION 
(FEET)1 

%  ANNUAL 
OCCURRENCE 
OF  SURGE 

DATE 

14.41 

5.0 

5.6 

November  14,  1981 

14.07 

4.4 

18.9 

February  4,  2006 

13.87 

- 

- 

December  1 1 ,  1973 

13.61 

5.1 

4.5 

March  3,  1999 

13.41 

5.3 

3.0 

December  3,  2007 

13.36 

- 

- 

December  3,  1982 

13.23 

3.8 

54.3 

December  1 , 2001 

13.21 

3.3 

>100 

January  2,  2003 

13.16 

- 

- 

January  27,  1983 

13.09 

- 

- 

February  7,  1978 

12.97 

- 

- 

January  18,  1973 

12.95 

2.8 

>100 

January  29,  1999 

1  dash  indicates  years  prior  to  hourly  water  elevation  data 
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Figure  2.  Peak  storm  surge  hydrographs  recorded  at  Toke  Point  Station  9440910. 


METHODOLOGY 

In  order  to  determine  flooding  risk  to  the  Shoalwater  Indian  Reservation  infrastructure  a  statistical 
analysis  was  performed  on  the  observed  water  levels  and  wave  heights  to  generate  the  boundary 
conditions  for  the  inundation  model. 


Statistical  analysis 

Extreme  water  levels  are  a  combination  of  astronomical  tide,  storm  surge,  and  wave  setup.  The 
tide  component  can  be  considered  a  deterministic  event.  Storm  surge  is  a  stochastic  event  and  must  be 
analyzed  statistically  using  observed  data  to  determine  the  frequency  of  occurrence  of  various  surge 
elevations.  In  this  analysis,  storm  surge  is  simply  defined  as  the  difference  between  the  predicted 
astronomical  tide  and  the  observed  tide.  Figure  3  shows  the  joint  distribution  of  tides  and  storm  surge 
from  1980  to  2007.  The  coincidence  of  high  tides  and  extreme  storm  surges  is  evident  by  data 
populated  in  the  upper  right  comer  of  the  figure.  The  frequency  of  high  water  levels  can  be  computed 
from  the  joint  distribution  using  the  convolution  method  originally  described  by  Pugh  and  Vassie 
(1978),  assuming  the  surge  and  tides  are  statistically  independent  of  one  another.  At  Grayeyard  spit  it 
is  well  documented  that  surface  gravity  waves  are  tidally  modulated  suggesting  storm  surge  is  also 
dependent  on  tide  and  the  convolution  method  may  not  be  applicable.  Instead,  a  simplified  approach  is 
utilized  to  develop  the  extreme  water  levels  simulated  in  the  inundation  model.  Representative  tidal 
elevations  of  mean  higher  high  water  (MHHW)  and  maximum  astronomical  tide  (MAT)  are  added  to 
the  storm  surge  extreme  distribution  to  specify  the  initial  water  level  conditions.  The  assumption  is 
wave  setup  is  minimal  at  the  tide  station.  This  is  reasonable  as  the  tide  gage  is  located  in  the  sheltered 
area  on  the  bayside  of  Tokeland  Peninsula.  An  extreme  storm  surge  occurring  at  MHHW  represents  a 
likely  occurrence  while  a  peak  storm  surge  coincident  with  MAT  is  less  likely,  but  represents  an  upper 
limit  to  the  extreme  water  levels  in  the  area. 

A  peak  over  threshold  analysis  is  performed  to  develop  the  storm  surge  extreme  distribution. 
Storm  surge  elevations  above  a  threshold  of  3.3  ft  yielded  NT  =  31  events,  or  approximately  one 
extreme  event  per  year.  The  extreme  distribution  following  a  Fisher-Tippet  type  I  (or  Gumbel 
distribution)  is  fit  to  the  data  using  the  least  squares  method.  The  form  of  the  cumulative  distribution 
follows 


F(x )  =  exp 


f 

-exp 

v 


x-B\ 

~4~J 


:  <  x  < 


(1) 
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where  parameter  A  governs  the  linear  scale  of  x  and  parameter  B  fixes  the  location  of  the  x-axis.  The 
distribution  is  fit  to  the  data  using  the  unbiased  plotting  position  formula  outlined  in  (Goda  2000). 
Figure  4  shows  the  fitted  distribution  with  the  90%  confidence  intervals.  The  storm  surge  return 
interval  (in  years)  is  calculated  following  the  relation 

R  =  A[1-F(xh)Y  (2) 

where  X  is  the  ratio  of  extreme  samples  to  the  number  of  years  of  data  (e.g.  X  =  NT/N),  and  F  (xu)  is  the 
probability  that  the  extreme  variate  x  does  not  exceed  a  given  value  xu  in  one  year. 

Numerical  wave  modeling 

Three  historical  cases  are  modeled  representing  the  three  highest  storm  surge  events  on  record 
since  1980.  Additionally,  six  hypothetical  storms  are  modeled.  Table  2  lists  the  nine  storm  cases 
simulated.  The  hypothetical  storms  represent  the  50,  2,  and  1%  annual  storm  surge  occurrence  (i.e.  the 

2,  50,  and  100  year  return  interval)  occurring  at  MHHW  and  MAT.  Extreme  storm  surge  events  are 
well  correlated  with  wave  height,  period,  and  direction  as  these  extreme  events  are  almost  always 
associated  with  south- southwesterly  tracks.  Extreme  water  levels  combined  with  large  offshore  waves 
result  in  flooding  to  the  Shoalwater  Reservation  shoreline.  Documented  damage  to  local 
infrastructure  has  typically  been  associated  with  dynamic  flooding  due  to  wave  runup,  overtopping, 
and  overland  wave  propagation  during  these  events. 

CMS -Wave,  is  a  two-dimensional  (2D)  spectral  wave  model  formulated  based  on  the  wave  action 
balance  equation  with  energy  dissipation  and  diffraction  terms  (Mase  2001;  Mase  et  al.  2005).  It 
simulates  a  steady-state  spectral  transformation  of  directional  random  waves  co-existing  with  ambient 
currents  in  the  coastal  zone.  The  model  operates  on  a  coastal  half-plane,  implying  waves  can  propagate 
only  from  the  seaward  boundary  toward  shore.  It  includes  features  such  as  wave  generation,  wave 
reflection,  and  bottom  frictional  dissipation.  The  model  utilizes  analytical  relations  to  include  wave 
setup,  runup,  and  overtopping  to  compute  inundation  depths  over  low  lying  backshore  regions  (Lin  et 
al.  2008). 

Data  Sources 

The  digital  elevation  model  used  to  construct  the  numerical  grid  is  generated  from  three  sources. 
Topographic  data  from  a  photogrammetric  survey  in  July  2008  are  used  to  generate  the  Graveyard  Spit 
condition  and  Tokeland  Peninsula  shoreline  and  backshore  regions.  Annual  Corps  hydrosurvey 
condition  survey  data  from  September  2008  of  the  north  entrance  channel  and  nearshore  area  are 
utilized  for  the  surrounding  bathymetry.  Finally,  the  (Taylor  et  al.  2008)  combined  topography  and 
bathymetry  tsunami  inundation  digital  elevation  model  is  utilized  for  coverage  outside  these  2008  data. 
The  model  grid  is  constructed  relative  to  the  mean  tide  level  (MTL)  tidal  datum. 

The  wave  modeling  is  performed  using  a  nested  grid  approach.  A  coarser  regional  scale  grid  with 
200  ft  cell  spacing  is  run  to  output  wave  information  to  “nesting  cells”  located  at  the  offshore  boundary 
of  the  higher  resolution  local  scale  grid.  This  40  ft  cell  spacing  grid  covers  Tokeland  Peninsula  and  the 
immediate  offshore  area.  The  two  grids  used  in  the  simulations  are  shown  in  Figure  5. 

Wave  spectra  for  the  cases  listed  in  Table  2  are  collected  from  the  Grays  Harbor  Coastal  Data 
Information  Program  (CDIP)  036  buoy  located  13  miles  northwest  of  the  Entrance  are  utilized  to 
specify  the  offshore  wave  boundary  condition.  For  the  case  of  the  March  3,  1999  storm,  there  is  a  gap 
in  the  CDIP  buoy  data;  therefore  the  spectra  from  the  National  Data  Buoy  Center  (NDBC)  buoy  46029 
offshore  the  Mouth  of  the  Columbia  River  located  43  miles  southwest  is  utilized.  In  this  case,  the 
wave  spectrum  is  transformed  using  the  maximum  entropy  method  to  the  130  ft  depth,  which 
corresponds  to  the  water  depth  at  the  offshore  boundary  of  the  regional  scale  grid.  Cells  representing 
land  are  specified  in  the  *. struct  file  to  compute  wave  setup  and  runup. 

Model  results  and  validation 

Model  results  indicate  that  the  central  breach  through  Graveyard  Spit  allows  the  most  wave  energy 
to  enter  North  Cove.  Alongshore  wave  height  from  Stations  1-22  are  shown  in  Figure  5  for  the  March 

3,  1999  storm  event.  In  the  current  condition,  a  maximum  significant  wave  height  of  Hs  =  2.8  ft,  or 
0.85  m  is  predicted  at  Station  18. 
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Figure  3.  Toke  Point  Station  9440910  joint  distribution  of  tide  stage  (ft,  MLLW)  and  storm  surge  (ft)  from 
1980-2007.  Color  scale  represents  number  of  events. 


Figure  4.  Storm  surge  elevation  return  interval  in  feet  at  Toke  Point  Station  9440910 
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Table  2.  Model  cases  run  in  CMS-WAVE 

Storm  Description 

Peak  Surge 
elevation 
(«) 

Total  water 
level 
(ft,  mllw) 

Wave 

height 

(«) 

Wave 
period  (s) 

Wave 

direction 

(deg) 

3-Mar-99  1000  UTC 

5.1 

13.6 

34.8 

16.67 

226 

3-Dec-07  1700  UTC 

5.3 

13.4 

31.5 

16.67 

214 

4-Feb-06  1300  UTC 

4.4 

13.9 

39.3 

15.38 

233 

2yr  @  MHHW 

3.9 

12.8 

26 

16.67 

226 

50yr  @  MHHW 

5.5 

13.4 

31.5 

16.67 

214 

100yr@  MHHW 

5.8 

13.9 

39.3 

15.38 

233 

2yr  @  MAT 

3.9 

14.9 

26 

16.67 

226 

50yr  @  MAT 

5.5 

16.5 

31.5 

16.67 

214 

lOOyr  @  MAT 

5.8 

16.8 

39.3 

15.38 

233 

1  MHHW  =  8.94  ft  MLLW 

2  MAT  =  11  ft  MLLW 

The  increase  in  wave  height  through  the  breach  brings  potential  adverse  impacts  such  as  larger  wave 
runup  on  the  shoreline  and  potential  for  greater  storm  related  flooding.  Figure  6  shows  ground  photos 
taken  by  members  of  the  Shoalwater  Reservation  and  demonstrates  flood  inundation  and  the 
consequences  of  overland  wave  propagation  and  debris  damage  associated  with  the  March  3,  1999  and 
February  4,  2006  storms. 

Inundation  Model 

Inundation  depth  is  defined  as  the  water  depth  associated  with  overland  flooding.  The  maximum 
water  depth  can  be  a  result  of  wave  runup  or  inundation  by  rise  in  water  levels  due  to  tide  and  storm 
surge.  The  total  wave  runup  consists  of  two  components:  (a)  rise  of  the  mean  water  level  by  wave 
breaking  at  the  shore,  known  as  the  wave  setup,  and  (b)  swash  of  incident  waves.  Wave  setup  is 
computed  as: 


dr/ 

dx 


pgh 


AS'  AS’ 


dx  dy 


dr/ 

dy 


pgh 


AS'  AS',, 


dx  dy 


(3) 


where  p  is  the  water  density  and  Sxx ,  Sxy,  and  Syy  are  radiation  stress  components  from  the  excess 
momentum  flux  caused  by  waves.  The  swash  oscillation  of  incident  waves  on  the  shoreline  is  a  random 
process.  A  statistical  parameter  describing  the  2%  exceedance  of  all  runup  is  defined  by  (Komar  1998) 
as: 

R2  =  2r|max  (4) 

In  CMS -Wave,  R2  is  calculated  at  the  land- water  interface  and  averaged  with  the  local  depth  to 
determine  if  the  water  can  flood  the  proceeding  dry  cell.  If  the  wave  runup  level  is  higher  than  the 
adjacent  land  cell  elevation,  CMS-Wave  will  flood  the  dry  cells  and  simulate  wave  overtopping  and 
overwash  at  them.  In  areas  which  have  been  flooded  to  allow  overland  wave  propagation,  the 
maximum  water  level  is  computed  as  the  wave  amplitude  plus  the  wave  setup,  or,  rj  +  0.5H1/3,  where 
H1/3  is  wave  height.  In  low  lying  regions  removed  from  overland  wave  propagation,  water  may  flood 
strictly  due  to  static  water  level  increases  by  storm  surge  and  tide.  Thus,  over  the  entire  model  domain 
the  maximum  flooding  depth  is  computed  using  the  following  relation: 
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Figure  5.  CMS-WAVE  model  domain  (a)  Coarse  grid;  (b)  nested  grid  with  observation  cell  locations  (water 
depth  in  meters,  MTL);  (c)  calculated  wave  height  in  meters  for  3  March  1999  10:00  UTC  simulation  at 
observation  stations. 
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Maximum  flooding  depth  = 

max  ( R2 ,  r|  +  0.5H1/3)  +  tide  elevation  +  surge  elevation  -  land  elevation  (5) 

In  shallow  water,  the  maximum  velocity  of  the  wave  (or  celerity)  is  only  a  function  of  water  depth  thus 
the  velocity  of  overland  wave  propagation  is  approximated  as 

V  =  Jgh  (6) 

where  g  is  gravitational  acceleration  and  h  is  the  local  water  depth.  This  is  a  simple  approximation  of 
the  velocity  as  the  true  velocity  will  be  affected  by  substrate  type  and  vegetation.  However,  the  intent 
of  this  analysis  is  to  develop  a  simple  metric  used  to  objectively  identify  potential  consequences 
associated  with  overland  wave  propagation  to  susceptible  infrastructure. 


Figure  6.  Inundation  of  Shoalwater  Indian  Reservation  uplands  and  damage  to  structures  (a)  3  March  1999 
storm  (b)  4  February  2006  storm  (c)  4  February  2006  storm 

PROPOSED  PLAN 


Dune  sizing 

The  dune  restoration  alternative  is  intended  to  restore  the  eroded  barrier  island  dune  system  with 
sand  dredged  from  the  adjacent  Willapa  Bay  north  entrance  channel.  The  crest  elevation  is  designed  to 
minimize  the  risk  of  wave  overtopping  and  maintain  a  sufficient  sand  reservoir  following  a  100  year 
storm  event.  In  the  narrowest  section  the  cross-sectional  area  of  the  dune  above  mean  high  water 
(MHW)  is  1785  ft2  per  lineal  foot.  The  restored  dune  is  12,500-feet-long,  with  a  top  elevation  of  +25 
feet  MLLW,  a  crest  width  of  20  feet,  and  a  side  slope  of  IV  on  5H.  In  the  updrift  reach  and  region 
where  the  breach  will  be  closed  a  protective  beach  berm  will  be  constructed  at  elevation  +15  ft 
MLLW.  Design  details  are  included  in  (US ACE  2009).  The  initial  dune  restoration  requires 
approximately  700,000  cubic  yards  (CY)  of  sand  dredged  from  the  entrance  to  Willapa  Bay.  It  is 
estimated  the  dune  will  need  to  be  renourished  every  five  years  with  approximately  250,000  CY  of 
dredge  material.  To  reduce  aeolian  transport,  the  leeside  crest  and  sideslope  of  the  constructed  dune 
will  be  planted  with  native  dune  grass. 
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Borrow  pit  sizing 

The  relative  water  depth  kh  of  the  incident  wave  conditions,  borrow  pit  cross  shore  width,  and 
borrow  pit  depth  have  the  greatest  potential  to  modify  the  surrounding  wave  environment  through  the 
processes  of  wave  reflection,  refraction  and  diffraction  (Michalsen  2004;  Michalsen  and  Haller  2008). 
To  minimize  impacts  to  wave  and  currents  the  borrow  pit  orientation  is  designed  to  follow  the  contours 
of  the  North  Entrance  channel  sideslope.  The  dredge  cut  will  be  approximately  8  ft  below  the  existing 
seabed  and  have  a  footprint  of  180  acres.  Since  2000  this  area  has  accreted  approximately  0.75  MCY 
per  year  as  the  channel  has  migrated  southward.  The  north  entrance  channel  presently  exhibits  depths 
of  90  ft  in  the  thalweg,  and  the  effect  of  the  borrow  pit  essentially  widens  the  channel.  Figure  7 
displays  the  effect  of  the  borrow  pit  and  has  relatively  minor  effects  on  the  wave  climate  (<  0.3  ft,  or 
0.1  m  change)  from  the  existing  configuration. 


RISK  ANALYSIS 


Without  project  condition 

Damage  to  structures  from  coastal  flooding  is  often  related  to  flood  depth  or  velocity  of  flood 
waters  (FEMA  2005).  High  velocity  flow  associated  with  wave  runup,  overtopping,  and  overland 
wave  propagation  has  the  potential  of  transporting  damaging  debris.  Potential  damages  associated  with 
flooding  can  occur  to  SR  105,  private  residences,  and  property  owned  by  the  Shoalwater  Indian 
Reservation.  A  structure  inventory  was  performed  for  all  property  located  within  the  Reservation 
Boundary  in  (US ACE  2009).  Flooding  depth  is  a  function  of  the  proximity  to  the  shoreline  and  the 
structures  first  floor  elevation  and  adjacent  grade.  The  threat  of  flooding  is  assigned  using  the 
following  classification  risk  matrix  presented  in  Table  3.  A  flooding  threat  to  a  structure  is  defined  as 
any  classification  exceeding  a  low  threat  (L).  As  a  validation  set,  the  March  3,  1999  storm  resulted  in 
54%  of  the  structures  classified  with  a  flooding  threat.  Figure  8  shows  a  detailed  inundation  map  for 
the  March  3,  1999  storm.  Documented  damages  by  local  residents  during  the  storm  confirm  model 
results  to  be  reasonable  at  predicting  inundation  depths. 

The  same  approach  is  utilized  to  describe  the  threat  during  the  50,  2,  and  1%  annual  storm  surge 
occurrence  intervals.  Table  4  lists  the  flooding  risk  of  Shoalwater  Reservation  structures  for  each 
occurrence  interval  at  MHHW  and  MAT,  respectively.  For  a  2%  annual  storm  surge  occurring  during 
MAT,  92%  of  the  structures  inventoried  are  classified  with  a  medium  or  high  flooding  risk. 
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Table  3.  Risk  assessment  tool  for  determining  consequences  to  inventoried 
infrastructure 


Risk  Matrix 


0 

a 

U) 

c 

'•u 

o 

o 


Q. 

0 

o 


D  >1.0 


.03<  D  <  1.0 


D  <.03 


M 


M 


V<1.0  1.0<V<5.7  V  >  5.7 


Velocity  of  Flow  (feet/second) 


Significance  of  Risk 


Medium  (M) 


Depth  of  flooding  in  structures  above  one  foot  causes 
significant  damage  to  structure  and  contents  and/or 
debris  damage  from  high  velocity  flows 


Flood  depths  inside  structures  up  to  one  foot  can 
cause  damage  to  structure  or  contents  and/or  debris 
damage  from  high  velocity  flows 


Low  (L) 


Nuisance  flooding  with  no  flooding  in  structures,  with 
still  water  that  poses  no  life  safety  threat.  No 
possibility  of  high  velocity  flow  or  debris  damage. 


1  velocity  is  a  function  of  depth;  therefore  high  velocity  flows  can  only  occur  with 
sufficient  depth. 


Table  4.  Percentage  of  Structures  at  Flooding  Risk  for  existing  condition  (Without  project)  alternative 
for  storm  surge  events  occurring  at  MHHW  and  MAT 

MHHW 

MAT 

Flood  and 
Storm 
Damage 
Risk 

50%  Annual 
Occurrence 

2%  Annual 
Occurrence 

1%  Annual 
Occurrence 

50%  Annual 
Occurrence 

2%  Annual 
Occurrence 

1%  Annual 
Occurrence 

Low 

85 

17 

12 

9 

8 

8 

Medium 

7 

40 

35 

35 

12 

12 

High 

8 

43 

53 

56 

80 

80 

With  project  condition 

North  Cove  will  continue  to  be  hydraulically  connected  to  Willapa  Bay  through  two  inlets. 
However,  the  restored  dune  will  significantly  reduce  the  frequency  and  magnitude  of  dynamic  water 
level  flooding  due  to  wave  runup  and  overland  wave  propagation.  Figure  9  depicts  the  new  risk 
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assessment  following  dune  restoration.  The  risk  to  the  inventoried  structures  is  reduced  from  54% 
(without  project)  to  7%  (with  project)  during  the  March  3,  1999  storm  condition.  Additionally  for  the 
6  hypothetical  cases,  flood  risk  is  significantly  reduced  (>70%  for  MHHW  and  >24%  for  MAT)  as 
displayed  in  Table  5.  During  storm  conditions,  the  restored  dune  will  mitigate  structure  damage  from 
debris  carried  inland  by  high  velocity  sheet  flows.  However,  instances  of  extreme  water  levels  caused 
by  large  storm  surges  at  MAT  will  still  flood  low  lying  topography,  thus  some  residual  flooding  risk  to 
structures  will  exist  even  with  the  restored  dune. 


Table  5.  Percentage  of  Structures  at  Flooding  Risk  with  restored  barrier  Island  condition  (With  project) 
for  storm  surge  events  occurring  at  MHHW  and  MAT 

MHHW 

MAT 

Flood  and 
Storm 
Damage 
Risk 

50%  Annual 
Occurrence 

2%  Annual 
Occurrence 

1%  Annual 
Occurrence 

50%  Annual 
Occurrence 

2%  Annual 
Occurrence 

1%  Annual 
Occurrence 

Low 

100 

87 

83 

83 

37 

32 

Medium 

0 

8 

6 

6 

19 

15 

High 

0 

5 

11 

11 

44 

53 

%  reduced 
to  Low 

15 

70 

71 

74 

29 

24 

CONCLUSIONS 

The  U.S.  Army  Corps  of  Engineers  has  determined  barrier  island  dune  restoration  as  the  most 
feasible  alternative  for  storm  damage  reduction  at  the  Shoalwater  Indian  Reservation  in  Willapa  Bay, 
Washington  USA.  Extreme  water  levels  coincident  with  strong  winter  storms  are  capable  of 
inundating  this  low  lying  topography  following  the  erosion  and  overwash  of  the  protective  barrier 
islands  known  as  Graveyard  Spit.  The  proposed  plan  would  dredge  material  from  the  Willapa  Bay 
north  entrance  channel  to  close  an  existing  breach  and  reconstruct  the  dune  on  Graveyard  Spit  to 
mitigate  flood  damages  to  the  Reservation  shoreline  and  uplands. 

To  quantify  the  benefits  of  the  project,  a  simple  tool  is  presented  for  identifying  flood  risk  to  the 
Shoalwater  Indian  Reservation  infrastructure.  Statistical  analysis  of  extreme  water  levels  and 
numerical  modeling  is  utilized  to  determine  the  extent  of  inundation  in  the  Reservation.  In  the  current 
condition,  it  was  determined  54%  of  the  inventoried  infrastructure  on  the  Shoalwater  Indian 
Reservation  is  at  risk  during  the  March  3,  1999  storm  event.  With  the  barrier  island  restoration  in 
place  it  was  found  that  this  risk  is  reduced  to  7%. 


12 


COASTAL  ENGINEERING  2010 


Figure  8.  Without  project  risk  assessment  for  Shoalwater  Reservation  infrastructure  (March  3,  1999  storm) 


Figure  9.  With  project  risk  assessment  for  Shoalwater  Reservation  infrastructure  (March  3, 1999  storm) 
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